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ABSTRACT 

In this paper, we consider the primary user detection problem in cognitive radio systems by using multi 
antenna at the cognitive radio receiver. An optimal square law combiner multi antenna based spectrum 
sensing technique is proposed using the multitaper spectrum estimation method. The multitaper spectrum 
estimation method (MTM) produces single spectrum estimate with minimum spectral leakage and variance 
using an orthonormal family of tapers, the Discrete Prolate Slepian Sequences (DPSS). An energy detector 
(ED), which is simple, but has a poor performance at a low signal to noise ratio (SNR). 

Keywords: Cognitive radio, spectrum sensing,multi taper spectrum estimation method, multi antenna 
spectrum sensing. ror. 



I. INTRODUCTION 

Radio Spectrum refers to the existing, natural 
medium that is used in different wireless 
communication systems and services: mobile, fixed, 
satellite-based, and low-power device 
communications systems (ultra wideband, sensors 
etc.)[l].Due to increase in wireless technologies , there 
is a problem of spectrum scarity. Cognitive radio (CR) 
is a new wireless communication technology that 
assigns a spectrum dynamically to the secondary 
users when the primary users are not using licensed 
spectrum. 

There are different types of spectrum sensing 
techniques are there like Matched filter, Energy 
detector ,cyclostationary ,multi taper 
method. Matched filter [2] and cyclostationary are 
classified as high formance spectrum sensing 
techniques but they require every primary user 
transmitted signal i.e which requires primary user 
information. The disadvantage in the cyclostationary 
detection [3] is which requires cyclic frequencies of 
the primary user and also it takes long sensing time. 
Similarly the practical implementation of Energy 
detector is easy but gives very poor performance at 
low signal to noise ratio (SNR). 



II. IMPORTANCE OF MTM METHOD FOR 
APECTRUM SENSING IN COGNITIVE RADIO 
NETWORK 

Multitaper spectrum estimation method (MTM) was 
proposed in 1982 by Thomson [4]. MTM uses an 
optimal bank of band pass filters (known as tapers or 
windows). These orthonormal tapers are called 
Discrete Prolate Slepian Sequences (DPSS) [6]. MTM 
produces a single spectrum estimate with minimum 
spectral leakage and good variance. The spectrum 
estimation in MTM is an approximation of the 
optimal estimate; the maximum likelihood (ML), One 
advantage for MTM compared to ML is the fact that 
it has lower computation complexity. Since the first 
development of MTM in 1982, this advanced method 
has been widely used in many applications. In 
addition to the power spectrum estimation in signal 
processing and communications applications, MTM 
is used in neurosciences, geophysics and sonar. 

III. THE IMPORTANCE OF USAGE OF MULTI 
ANTENNA BASED SPECTRUM SENSING IN CR 
SYSTEMS 

In order to improve the spatial diversity, classical 
wireless communications use multi antenna at 
transmitter (Tx) or receiver (Rx) or both. Such 
diversity improvement increases the system data rate 
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and capacity. The reason behind this is that as the 
distance between antennas is chosen properly there 
will be a high probability of receiving independent 
fading through these different antennas [9], 
Therefore, the fading effects will be mitigated. 

IV. ENERGY DETECTOR 

The probability of detection P d is defined as is the 

probability that the CR detector decides correctly the 
presence of the PR's signal. Therefore, this type of 
probability is related to the P(Z);//j) which 
represents the PR's signal plus noise 

00 

P d =P r {D>y/H l ] =\p{D;H x )dt (l) 

r 

where P r represents the probability .It is clear 

that p d is the integration over f(^D; H ,) from 
the threshold y to infinity. Since the 

P\D\ H ^ ) is assumed as Gaussian distribution; 

Pd can be defined finally as follows[5] 

y-E[p(D;/ H Y ) 



(2) 



yvAflP(Affi)] 

where E and Var, are the mean and variance of the 
given distribution respectively. The term <2(<r) 1S 
the complementary cumulative distribution function, 

.2 

-a GO t 



dt 



(3) 



The probability of false alarm, P f is defined as the 

probability that the CR detector decides by mistake 
the presence of the PR"s signal .Therefore, this type 
of probability is related to the ,P{D\H 0 ) 
which represents the noise only. 



au 

P f = P r {D > y I H 0 } =\p(D;H 0 )dt 



(4) 



In this case, it is clear that the P* is the integration 

over from the threshold y to 

infinity[5] 

r-grtDi/ir.) 

jVAHP(D;H<,)] 

where E and Var, are the mean and variance of the 
given distribution respectively. 



V. MULTI TAPER SPECTRUM ESTIMATION 
METHOD 

some of the advantages of multi taper estimation 
method are 

1. MTM is an energy based spectrum sensing 
technique 

2. MTM is a wideband spectrum sensing technique 

3. MTM does not need any prior information about 
the PR"s signal (i.e., non coherent nor partial 
coherent). 

4. MTM needs to know the noise variance to control 
the threshold. 

5. MTM minimizes the spectral leakage outside the 
band and improves the variance of estimate. 

A general procedure for multi taper spectrum 
estimation consists of four steps. First, choose a time- 
bandwidth product CO = NW where N is the sample 
size for estimation, and W is the bandwidth 
normalized by the sample rate. The choice of 
parameter CO is a trade-off between spectral 
resolution and variance. Typically, 3 < CO < 6, and we 
use K = 2C0-1 data tapers to perform the estimation 
for the reason that will be given below. Second, 
compute the data tapers, i.e., Slepian sequences. 
Denote the kth data taper in vector form by wk. It can 
be calculated by the following Eigen equation 



KW k = W k 



(6) 



where the (i, j) th entry of Toeplitiz matrix K is 
defined by 

■ •- ' N ( 7 ) 



K 



i, 7 = 1,2,3. 



The Eigen values A K range between 0 and unity, 
and are organized in an descending order such 

that A l > A 2 > A N . The first 

K ~ [2NW] of them are dominated, close to unity, 
whereas the rest are negligible. Moreover, the tapers 
of lower order have much stronger energy- 
concentration capability than their high-order 
counterparts, suggesting that it suffices to use the 
first k tapers for spectral estimation. Third, the 
corresponding eigenspectra are defined by Fourier 
transform of the tapered (i.e., windowed) data 
sequences, resulting in 

N 

Y k (f) = Y J w k {n)x{n)e- i2Kf ' 1 k = 0,1,2 it -1(8) 

as a function of frequency f, Here, wk(n) signifies the 
nth entry of wk. Given a finite sample-size constraint, 
the energy distribution of each eigenspectrum falls 
between the bandwidth from f -W to f +W, due to 



LIRA I 2014 I Volume 1 I Issue 4 



Pag 



164 



International Journal of Research and Applications (Oct - Dec © 2014 Transactions) 



International Conference on Emerging Trends in Electronics & Telecommunications (ICETET-14) 



energy-concentration property of Slepian sequences 
Finally, we combine the Eigen spectra to form a 
single spectrum estimate, 

* A (/)=!>* I W)f (9) 

k=0 

where the weighting factor 

0t k = X k I(Xq +A l + + is used to account for 

the relative importance of the Eigen spectrum of 
concern. 

The received PR signal, at the CR receiver, is 
sampled to generate a finite discrete time samples 
series {x lm ; t = 0,1,2, N-l, m = l,2, M 

where m denotes the antenna number and t is the 
time index The discrete time samples are dot 
multiplied with different tapers v tk (N,W) (Tapers 

are DPSS). 

The associated Eigen value of the k th taper 
is A k (N,W) . The product is applied to a Fourier 

transform to compute the energy concentrated in the 
bandwidth (-W, W) centred at frequency, f. The half 
time bandwidth product is NW. The total number of 
generated tapers is 2NW. 

The binary hypothesis test for CR spectrum sensing 
at the l' h time, and using the m' h antenna branch is 
given by 

H l :x lim (l) = s t (D + w, m (l) (10) 

where 1=0, 1, 2 L-l OFDM block's index, and 

x t, m ( l ),w [m (l) and s t (l) denote the cr 

received, noise at the branch m and PR's transmitted 
samples. 

The transmitted PR signal is distorted by the zero 

mean AWGN, w t m (I) ~ N(0, a w 2 ) at the output 

from the different antenna branches, which are 
independent and with identical variance. For K 
orthonormal tapers used in the MTM, there will be K 
different Eigen spectrums produced from each 
antenna defined as [6] 

Y k Mi) = Xv.^WWXx,^/))^^* (11) 



where // ~ n ■> N ■> N 
normalized frequency bins. 
The power spectrum estimate given by [4] 



N-l 

N are the 



Y J ^ K {N,W)\Y KM {f, ) 
S m ( f ) = ^° 



^A k (N,W) 



(12) 



4=0 



The energy detector, when the samples are taken at 
uniform time spacing, gives the power spectrum 
density estimation[7] 

2 



'ED \J i 



(fi) = 



N-l 



t=0 



(13) 



A) Decision statics for the MTM and ED 

The decision statistic over L using MTM is defined 

for the TTl antenna as 



K-l 


N-l 






k=0 


r=0 



/=0 



(13) 



The decision statistic over L using ED is defined for 
the m th antenna as 

1=0 f=0 



(14) 



B) Mean and Variance of MTM and ED using Single 
Antenna 

For single antenna MTM-based spectrum sensing, 
and according to the central limit theorem, if the 
number of samples L, is large, the decision statistic, 
DEC'^j-j^ifj) is asymptotically normally 
distributed with mean[8] 

2 



EiDEC'^im 



\lK(E s +<j w 2 ) 



H 



0 



(15) 



and variance (VAR) 

VAR(DECl m (f i ))=l / \ 

[2LC 4<7 ff ((T v +2E S ) H t 

(16.1) 



c = E 



k=0 (17) 
The decision statistic, for the Energy Detector 
DEC ™ D (fj) is asymptotically normally distributed 
with mean[9] (18) 
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and variance (VAR) 



VARiDEC'M))- 



\2L(7 W 

\2La 2 (a 2 +2E S ) 



H 0 



(19) 

C) Detection and False Alarm probabilities 

For a normally distributed decision statistic, 

DEC {ff) the probabilities of detection P d (fj) and 

false alarm, P f (/. ) are defined as 

P d (f i ) = P(DEC(f i )>ylH l ) 

y - E(DEC(fj )/H 1 ) (20) 



^VAR{DEC{f i )/H 1 ) 



PJf i ) = P(DEC(f i )>y/H Q ) 



y-EiDECif^/H,) 



(21) 



\JVAR(DEC(f;)/ H 0 ) 
The probability of miss detection P mi (fj) is defined 



as 



P im (f,) = P(DEC(f l )<y/H 1 ) 

= l _ Q( _r-E(DEC(f i )/H l ) 



t]VAR(DE C(f t )/H 1 ) 
The term Q(^) is the complementary cumulative 

distribution function = [ e 2 dt and Y 

represents the chosen threshold. 

Thus, the different probabilities of MTM using single 

antenna can be redefined as [8] 

2, 



tMTM , r n. y-LK(E s +a„) 



P d (f i ) = Q( 



yj 2LC 1 A, 1 ,g w 2 (a w 2 + 2E S ) 



) 



E[DEC m ED {fi)] = 



\La w 

L(E s +a 2 ) 



(23) 



(24) 



y-LKa„ 



y}2LC 2 A z a w 4 ) 



(25) 



P,r(f,) = l-Q(- r LJ£B^ > =) (26) 
yJ2LC 2 A^(T w (cr w +2E S ) 

The number of samples required by MTM using 
single antenna ( E M ™ ) can be written as[8] 



r MTM 



f MTM j MTM \ 

a — b 



v KE s 
where a and b are 



(27) 



mtm Inr/^ 2 ! ^ A n~ x ( t> m ™ ( f \\ 
a =y]2LC A^(J W Q (P f (/,.)) 



and 



r , \MLKa 2 H n 

e[dec mtm _ slc (/,.)] = 

[MLK(E S +(T W 2 ) H 



b Mm = ^2LC 2 ^a w \cJ w 2 +2E s )Q-\P d Mm (/,)) 



E[DEC ED _ SLC (f i )\ = -\ 

[ML(E s +cr w 2 ) H l 

D) Mean and Variance of MTM and ED using M 
Antennas 

The decision statistics for square law combining 
using M antennas for both techniques MTM, and ED 
respectively as follow [12] 



DEC, 



MTM-SLC 



L-l 



^v, t (iV,H').v,,„(;)^ i2 * 



(22) and 



m=l ;=o 



M L-l N-l 



(28) 
(29) 



DEC ED _ SLC (f,) = Z j ^^ JDc ^' 

=1 1=0 t=0 

The decision statistic using square law combining is a 

sum of identical and independent normally 

distributed M antennas decision statistics. Thus 

DEC mtm-slc (fi^ the mean of the using M 

antennas [13] can be definedas (30) and the variance 

\2MLC 2 Lo' H n /o-n 

VAKDEC Mm _ s JM = \ 1 \ 0 ( 31 ) 

For the ED, the mean of the decision statistic, using 
square law combining through M antennas can be 
defined as (32) 

VI. SIMULATION RESULTS 

In our system each node of the cognitive radio (CR) 
network uses ,64-FFT with sampling frequency 
20MHz.The primary user'(PR) transmitter uses 64- 
IFFT with symbol duration Ts=0.05ps, and transmits 
QPSK signal with normalized equal to 1 with each 
subcarrier. The multi antenna technique MTM-SLC, 
ED-SLC. In MTM techniques the used half-time and 
width product is NW=4, and the number of tapers 
k=5.In fig. land fig 2 consider the AWGN channel 
with SNR=-10dB and 20 OFDM blocks (i.e L=20) 
used in sensing and the number of samples used is 
(L=20) X (N=64) =1280, which approximately 
corresponds to 64 ps. 
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Fig. 1. Probability of detection Vs probability of false 
alarm using MTM and ED with M=4 Antennas at the 
Receiver. 



Fig. 2. Signal to noise ratio (SNR) Vs probability of 
detection for M=4 Antennas at the Receiver. 




Fig. 3. Signal to Noise ratio (SNR) Vs Number of 
samples (dB) for M=4 Antennas at the Receiver. 
VII. CONCLUDING REMARKS 
Here we consider the detection of primary user by 
taking MTM and ED detectors. Multiple antennas 
(M=4) used at the receiver side. Using multi antenna 
in MTM-SLC gives more improvement in 
performance compared to that for ED-SLC and also 
MTM method requires less number of samples to 
detect the primary user compared to the ED. The OR 
rule cooperation is used to cooperate the generated 
binary decisions from the individuals CR nodes at a 



main CR-BS. Then the final decision is declared to 
the CR nodes. 
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Where A-^ is 

_ ^X k \N,W) + 2X 0 {N,W)X i (N,W) + U 0 (N,W)A 2 (N,W) + .+ 2X 0 {N,W)X KA {N,W) + 2\{N,W)X 2 (N,W) + 

1 h 2l x (N,W)^ {N,W) + .+ 2\ (N,W)A KA (N,W) + + 2X K _ 2 {N,W)X KA (N,W) 
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